Scaling behavior of the energy gap of spin-| AF-Heisenberg chain in both uniform and 

staggered fields 
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We have studied the energy gap of the ID AF-Heisenberg model in the presence of both uniform 
(H) and staggered (h) magnetic fields using the exact diagonalization technique. We have found 
that the opening of the gap in the presence of a staggered field scales with h v , where v = u(H) is 
the critical exponent and depends on the uniform field. With respect to the range of the staggered 
magnetic field, we have identified two regimes through which the if-dependence of the real critical 
exponent u(H) can be numerically calculated. Our numerical results are in good agreement with 
the results obtained by theoretical approaches. 

PACS numbers: 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 

The effect of external magnetic fields in the quan- 
tum properties of low-dimensional magnets has been of 
much interest in recent years. Experimental and theoret- 
ical studies of these systems have revealed a plethora of 
quantum flactuation phenomena, not usally observed in 
higher dimensions. The magnetization processes in an- 
tiferromagnetic (AF) spin chains and ladders have been 
under intensive investigation using novel numerical tech- 
niques. The progress in the experimental front is achivcd 
by introduction of high-field neutron scattering studies 
and synthesis of magnetic quasi-one dimensional systems 
such as the spin-i antiferromagnet Cu benzoatei^ and 
Yb^Asgr^^- Due to these developments we can now ob- 
serve the effect of a staggered magnetic field (or even 
more complicated interactions) on the low energy behav- 
ior of a one-dimentional quantum model in the labora- 
tory. 

There exist different mechanisms for generating a stag- 
gered field in a real magne1jZ»£»2. In Cu benzoate the al- 
ternating crystal axes is the source of such a field. Dender 
et.ali showed that an effective staggered field can be gen- 
erated by the alternating g-tensor. Theoreticaly, Afflec 
et.alii have studied how an effective staggered field is gen- 
erated by Dzyaloshinskii-Moriya (DM) interaction if the 
crystal symmetry is sufficiently low. They showed that 
in the presence of DM interaction along the AF chain, an 
applied uniform field T? generates an effective staggered 
field h . Ignoring small residual anisotropics, they ob- 
tained an effective hamiltonian where a one-dimensional 
Heisenberg AF chain is placed in perpendicular uniform 
(H) and staggered (h) fields 

H = Y^J^j^i+i - HS* + h{-iyS?] (1) 



It is expectec&iS that the staggered field induces an 
excitation gap in the S — Heisenberg antiferromagnetic 
(AF) chain, which should be otherwise gapless. Such as 
excitation gap caused by the staggered field is indeed 
found in real magnetsi&ii. 



In the absence of the staggered magnetic field (h = 0) 
and the uniform magnetic field (H = 0), the spectrum 
is gapless. In the ground state, the system is in the 
spin-fluid phase, where the decay of correlations fallow 
a power low. When a uniform magnetic field is ap- 
plied the spectrum of the system remains gapless until 
a critical field H c — 2 J, is reached. Here a phase tran- 
sition of the Pokrovsky-Talapov typei^ occurs and the 
ground state becomes a completely ordered ferromag- 
netic stated. Since the uniform magnetic field does not 
destroy the exact integrability of the Heisenberg model, 
the eigenspectra is exactly solvable. Applying a stag- 
gered magnetic field, the integrability is lost. The ap- 
plication of a staggered magnetic field when H = 0, 
produces an antiferromagnetically ordered (Neel order) 
ground-state and induces a gap in the spectrum of the 
model. Heisenberg model in both staggered and uni- 
form fields has been recently studied^ using density ma- 
trix renormalization group (DMRG). It is shown that 
bound midgap states generally exist in open boundary 
AF-Heisenberg chains. The gap and midgap energies in 
the thermodynamic limit are obtained by extrapolating 
numerical results of small chain sizes up to 200 sites. It 
is revealed that some of the gap and midgap energies for 
the half-integer spin chains fit well to a scaling function 
derived from the quantum Sine-Gordon model, but other 
low energy excitations do not fit equally well. 

In this paper, we present the numerical results ob- 
tained on the low-energy states of the ID AF-Heisenberg 
model in both uniform and staggered fields using an exact 
diagonalization technique for finite systems. We calculate 
the spin gap as a function of applied staggered field in 
the presence of small uniform field (0 < H < 0.1). With 
respect to the range of the staggered magnetic field, we 
show that there are two regimes in which we can com- 
pute the real critical exponent of the energy gap and 
it is important to note to which one of these regimes 
the numerical data are related. In Sec. II we discuss 
the scaling behavior of the gap using the available lim- 
iting behaviors. The leading exponent of the staggered 
field h, depends on H boath in finite size and thermo- 
dynamic limit. In Sec. Ill, we explain how, in certain 
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limits the numerical calculations may produce incorrect 
result for the critical exponent. We apply a perturbative 
approach 25 to find the correct critical exponent in the 
small- X (x = Nh v ( H ') regime. In Sec. IV, we increase 
the scaling parameter x and find the correct critical ex- 
ponent in the large- a; regime. Finally, the summary and 
discussion are presented in Sec. V. 



II. THE SCALING BEHAVIOR OF THE GAP 



In the high field neutron-scattering experiment on Cu 
benzoatei, which is a quasi-one dimensional S = s an- 
tiferromagnet, the magnetic field induces a gap in ex- 
citation spectrum of the magnet. The observed gap is 
proportional to Hq 65 , where Hq is the magnitude of 
the applied field. This exponent of about | describ- 
ing the field dependence of the gap obtained in differ- 
ent experiment a 6 ' identify the source of this gap as the 
staggered field. 

Using bosonization techniques, Affleck et.al showed 
that, the gap scales as 



A(h,H) ~ h u{H) 



(2) 



where v{H ) is the critical exponent of the gap and when 
H is stricly 0, v(H = 0) = |. The iJ-dependence of the 
exponent v{H) is studied numerically in Ref. [15]. Their 
approach is based on the 77-exponent, defined through the 
static structure factor of the model in the absence of a 
staggered field (h — 0). They show that there is a relation 
between the critical exponent of the gap and 77-exponent. 
Then by computing the 77-exponent of the structure fac- 
tor of the model, they predict the independence of v{H). 
Similarly, In an interesting recent worki£, the effect of an 
external field on the gap of the 2D AF Heisenberg model 
with DM interaction has been studied. It is shown that 
the effect of the external field on the gap can be predicted 
by investigating the on-site magnetization of the model. 

Here we study the evolution of the gap, using the con- 
formal estimates of the small perturbation h <C 1, and 
the finite size scaling estimates of the energy eigenvalues 
of the small chains in the presence of the staggered field 
(h ^ 0). We argue that there are two regimes in which 
the real critical exponent can be numerically calculated 
and it is very important to note to which one of these 
regimes the numerical data are related. 

Let us rewrite the Hamiltonian (1) in the form 

H = H + V 

Ho = ^[Jjj.jj+i — HSj] 



v = hj2(-iys>, 



(3) 



where Hq is exactly solvable by the Bethe ansatz and 
the staggered field ft < 1 is very small. For a small 
perturbation V, we can use conformal estimates. The 



large distance asymptotic of the correlation function of 
the model in the absence of the staggered field (h = 0) is 
obtained^ as 



-l) n 



(4) 



Where a{H) is a function of the uniform (H) field and 
is found using the Bethe ansatz as 



a(H) ~ 1 



1 



21n(l/i?)' 



H -> 



(5) 



where a(0) — 1 and a(2) = ^. By investigating the per- 
turbed action for the model and performing an infinites- 
imal renormalization group with a scale A, one can show 

that the staggered magnetic field scales as bl = h\ 2 ~ 
which leads that the energy gap scales as Eq.(f5|) by crit- 
ical exponent 



v(H) = 



a(H)' 



(6) 



which is also obtained with the bosonization technique in 
Ref. [7]. For example, in the absence of a uniform mag- 
netic field, A ~ ft,3, in agreement with the bosonization 
and experimental results. Since, by increasing the uni- 
form field H, a(H) decreases, thus we can conclude that 
the critical gap exponent v(H ) drops with increasing uni- 
form field H. 

To make a numerical check on effect of the uniform 
field on the energy gap we have implemented the mod- 
ified Lanczos algorithm 18 on finite-size chains (N = 
12, 14, 24) using periodic boundary conditions to cal- 
culate the energy gap. We have computed the energy 
gap for different chain lengths in the cases of the uniform 
fields < H < 0.1. The energy gap as a function of the 
chain length (N), uniform (H) and staggered (h) fields 
is defined as 



A(N,h,H) =Ex(N,h,H) -E (N,h,H), 



(7) 



where E is the ground state energy and E\ is the first 
excited state. In the absence of staggered field (h = 0) 
, the spectrum of the AF Heisenberg model is gapless 
up to H = 2 J. The gap vanishes in the thermodynamic 
limit proportional to the inverse of the chain length^ 



lim A(N,h = Q,H)- 



A(H) 



N- 



N 



(8) 



The coefficient A is known exactly from the Bethe 
ansatz solution 2 ^ and also can be computed in principle 
by the methods of conformal invariance and finite-size 
scalin g 21 i 22 ' 23 . 

When the staggered field is applied, a non zero gap de- 
velops. Thus the staggered field h c = is a critical point 
for our model. In general, the critical point of an infinite 
system is defined, in the Hamiltonian formulation, as the 
value of h, h c , at which the mass gap A(h, H) vanishes 
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as Eq.©. With our Lanczos scheme we can compute 
A(N,h,H), which approaches A(h, H) when N is large. 
The natural measure of the deviation of the finite sys- 
tem from the infinite one is -k-, where L is the linear 
dimension of the finite system (L = Na, a is the lattice 
spacing) and Lq is the correlation length of the infinite 
system (Lq — £a). Thus, we assume that A(N,h,H) 
depends on h through -k- as 



A(N,h,H)~f(^) = f(x) 



(9) 



where x — Nh u ^ H ^ is a scaling parameter, and f(x) is 
the scaling function. As expected, this equation behaves 
in the combined limit 



N 



oo, 







[x > 1) 



(10) 



as Eq.©, thus we assume the asymptotic form of the 
scaling function f(x) 



fix) ~ x* 



(11) 



In addition, we need a factor to cancel the N-dependence 
of f(x) as TV — ► oo. This factor must be of the form 
AT -1 . Thus, we have 

A(N, h, H) ~ Ar _1 /(x), (12) 

If we multiply both sides of Eq. (fT2|) by N we get 

lim NA(N,h,H)~x. (13) 

JV— >oo(x»l) 

Eq.lJTSJ) shows that the large- a; behavior of NA(N, h, H) 
is linear in x where the scaling exponent of the energy 
gap is v{H). 

III. SMALL-X REGIME 

Since the scaling of the gap can only be observed in the 
thermodynamic limit and for very small value of h, we 
compute the energy gap of the model for several values 
of staggered field 0.001 < h < 0.01 and different chain 
lengths N = 12, 14, 24 for fixed uniform fields < H < 
0.1. We have plotted the values of NA(N,h,H) versus 
Nh"W for H = 0,0.03,0.05,0.07,0.09. The results have 
been computed on a chain with the periodic boundary 
conditions. According to Eq. (|13p . we have found from 
our numerical results that the linear behavior is very well 
satisfied by v(H) = 2.0 independent of H. This is very 
far from the correct value of critical exponent v(H) < | 
(Eq.©). 

Note that the horizontal axes values in the small- x 
regime are limited to very small values of x = Nh v ^ H ^ < 
0.0024. Thus, we are not allowed to obtain the real scal- 
ing exponent of the gap which exists in the thermody- 
namic limit (N > oo or i > 1). 

When x is small, or in other words h is very small, 
we might be away from the thermodynamic behavior to 




FIG. 1: The value of scaling function gi(N,H) at the 
fixed uniform field H = 0, versus the chain length N = 
12, 14, ...,24. The best fit is obtained by v{H = 0) = 
2.05 ± 0.01. In the inset, the function gi(N,H) is plotted 
versus N at the uniform field H = 0.05 and the best fit is 
obtained by v(H = 0.05) = 2.33 ± 0.01. Data for different 
staggered fields 0.001 < h < 0.005 fall exactly on each other. 



observe the correct scaling. For very small h in the finite 
size systems (N ~ 24) the value of x will be small (x <C 
1), which avoid us to get the information on the large-x 
behavior of the scaling function f(x). In this case, the 
values of the energy gap coming from a finite system are 
basically representing the perturbative behavior 2 ^, which 
we reproduce for the convenience . 

We start from the Hamiltonian Eq.©. The energy 
eigenstates of the Hq carry momentum p = or p = tt 



T | ip n (h = 0,H)} = ±\ ip n (h = 0, H)), 



(14) 



where, T is translation operator and {| ip n (h = 0, H))} 
are eigenstates of unperturbed Hamiltonian Hq. Since 
the operator $^,-(— l)- 7 ^ changes the momentum of the 
state by tt, we obtain that 



(V>«(0,#)| V| ^n(0,if))=0 



(15) 



Thus, the gap can be rewritten in the following form 

A(N,h,H) = A{N,0,H)+ 9l (N,H)h 2 + 
••• +g n (N,H)h 2n , 



(16) 



where n is an integer. It is a good approximate to neglect 
the effect of higher order terms for h < 0.01. Because 
the second-order perturbation correction is not zero in 
the staggered field, the leading nonzero term is h 2 . If 
the small-x behavior of the scaling function is defined as 
f(x) ~ , we find that 



v(H)<f>, 



(17) 
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This shows that in the small- x regime, NA(N, h, H) is a 

2 

linear function of x vl - H K This is in agreement with our 
data in small- a: regime, where <fr s = 1 and according to 
Eq.([l7]), the value of v{H) is found to be u(H) = 2.0. 

Let us consider the large- TV behavior of g\ (N, H) at 
fixed H as 



lim gx{N,H) 

iV— >OC 



(18) 



This leads to 



A(N, h, H) 



A{H) 
N 



(l + b 1 (H)N^ H '>+ 1 h 2 ), (19) 



where b\{H) is a constant (at fixed H). We can write 
Eq. (p~9|) in terms of the scaling variable x 



NA(N, h, H) ~ A(H)(1 + N^ H)+1 ~^m x ^m ) (20) 

For large- N limit this equation should be independent of 
N, leading to the relation between ji(H) and v{H) as 



u{H) = 



1 



(21) 



The above arguments propose to look for the large-TV 
behavior of gi(N,H). To determine the n(H) exponent, 
we have plotted in Fig.l the following expression versus 
N 



9i(N,H) 



A(N,h,H)-A(N,0,H) 
h 2 



(22) 



for fixed values of staggered field h (0.001 < h < 0.005), 
and different sizes, N — 12, 14, 24 at the uniform field 
H = 0. We found the best fit to our data for fj,(H = 0) = 
2.04±0.01. The inset in Fig.l shows the gi(N,H) versus 
N at fixed H = 0.05. In this case, the best fit, found 
for (i(H = 0.5) = 2.33 ± 0.01. Our data for different 
h values, fall perfectly on each other, which shows that 
our results for gi(N,H) in fixed uniform field H, are 
independent of the staggered field h as expected. By 
using Ea.pi]) we have found, v{H = 0) = 0.66 ±0.01 and 
v{H = 0.05) = 0.60 ± 0.01. We have also implemented 
our numerical tool to calculate the critical exponent i>{H ) 
at H — 0.03,0.07,0.09. The results have been presented 
in Table I. 



IV. LARGE-X REGIME 




FIG. 2: Difference between the two lowest energy levels and 
the ground state energy as a function of the staggered mag- 
netic field for finite chain length N = 24 and H = in the 
region 0.01 < h < 0.4. 



< 
Z 



N=18 
N=20 
N=22 
N=24 




Nh" 



FIG. 3: The product of energy gap and chain length (ATA) 
versus Nh u{H) at the uniform field H = 0. In the range of 
staggered field 0.04 < h < 0.4, linear behavior is obtained 
by choosing v(H = 0) = 0.66 for all different chain lengths 
N = 18,20,22,24. In the inset, iVA is plotted for uniform 
field H = 0.05. The linear behavior is obtained by v(H = 
0.05) = 0.63. Data for different chain lengths fall on each 
other. 



In our numerical calculations because of memory is- 
sues, we were limited to consider the maximum chain 
length N = 24. Therefore the value of x cannot be 
increased by increasing the size of chain. The prob- 
lem appears if the calculation is done by density ma- 
trix renormalization group (DMRG) method. In that 
case, we may extend the calculation to larger sizes, 
N ~ 200, which cannot increase x much larger than one 



(for 0.001 < h < 0.01). On the other hand we may in- 
crease the staggered field for increasing x. But we should 
note that, in general there is usually level crossing be- 
tween the energy levels in finite size systems. Which 
can change the behavior of the gapSi and lead to incom- 
mensurate effects. As an example, the dependences on 
the excitation energies of the three lowest levels on the 
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FIG. 4: A graph of the critical exponent v versus uniform 
field H . Both critical exponents v a (squares) and vi (circles) 
start at the known value | and then drop with the increasing 
uniform field H. In the inset, we have plotted the critical 
exponent v which is obtained directly by extrapolating the 
numerical results of the energy gap in the range of the stag- 
gered field 0.001 < h < O.Of. 



TABLE I: The exponent of gi(N, H) and gap exponent versus 
different values of the uniform field H in small-x and large-x 
regimes. 



H 



VL 



0.0 
0.03 
0.05 
0.07 
0.09 



2.05 
2.20 
2.33 
2.58 
2.80 



0.65 
0.63 
0.60 
0.56 
0.53 



0.66 
0.64 
0.63 
0.62 
0.60 



v the plot will not be linear in x. Also, our data for 
different size chains fall perfectly on each other, which is 
expected from the scaling behavior. 

We have extended our numerical computations to con- 
sider values of the uniform field in the small region 
< H < 0.1. The results have been presented in Ta- 
ble I. We have listed /i, the resulting v s that is obtained 
from perturbative approach, and result of the large-x 
regime vi for different values of the uniform field H. In 
Fig. 4 we have plotted the critical exponent v(H) ver- 
sus the uniform field H. As it is clearly seen from this 
figure, v s and vi start at the known value | and then 
drop with the increasing uniform field H . The expo- 
nents are in good agreement with each other and show 
good agreement with the exponents derived in the field 
theoretical approach (Eq.©). The inset in Fig. 4 shows 
the H-dependence of the critical exponent v which is ob- 
tained directly by extrapolating the numerical results 
of the energy gap in the range of the staggered field 
0.001 < h < 0.01 . It is clearly seen from this figure 
that the behavior of the gap in finite systems is different 
from its behavior at the thermodynamic limit, and with 
respect to the range of the staggered magnetic field the 
behavior of the gap deviates from the predicted scaling 
behavior. 

On the other hand, Fouet et.al studied^ the gap- 
induced by the staggered field h at the saturation uniform 
field H c = 2J. Using field theoretical arguments, they 
found that the gap scales as A(h, H c ) ~ ft 4 / 5 . By apply- 
ing the DMRG method for system sizes up to N = 100, 
they have also computed the exponent of the energy 
gap 0.81. We have extended our numerical computa- 
tions to consider the uniform field at the saturation value 
H c . In this case, we have obtained, v s (H c ) = 0.78 and 
vi(H c ) — 0.82 in good agreement with the Foet results. 



V. CONCLUSIONS 



staggered field h (0.01 < h < 0.4) are shown in Fig. 2 
for N — 24 and H = 0. From this figure, it can be 
seen that the two lowest excited states do not cross each 
other in this case. This means that we can increase the 
value of x by increasing h up to h = 0.4. Since that 
the regime where we can observe the right scaling of the 
gap is in very small value of h, we have restricted our 
numerical computations to upper limit of staggered field 
h = 0.4. In this respect, we have performed some numer- 
ical computations on the Hamiltonian Eq.{T]) for large- a; 
values, and the results are plotted in Fig. 3. In this fig- 
ure the range of the staggered field is 0.04 < h < 0.4, 
which causes that we get large- a: values for the consid- 
ered chain lengths N = 18, 20, 22, 24 and H = 0. The 
inset in Fig. 3 shows NA(N, h, H) versus x = Nh v{H) at 
fixed uniform field H = 0.05. In this case, we have ob- 
tained, v(H = 0) = 0.66 and v(H = 0.05) = 0.63. It 
should be mentioned that if we choose another value for 



In this paper, we have studied the energy gap of the 
ID AF-Heisenberg model in the presence of both uni- 
form (H) and staggered (h) magnetic fields using the 
exact diagonalization technique. We have implemented 
the modified Lanczos method to obtain the excited state 
energies with the same accuracy as the ground state one. 
We have been limited to a maximum of N = 24, because 
of memory considerations. We have shown, if the energy 
gap in the thermodynamic limit had been obtained by ex- 
trapolating the numerical results for finite systems, then 
the behavior of the gap may deviate from the predicted 
scaling behavior. This deviation depends on the range 
of the staggered magnetic field (h). We have found in 
the range of very small values of the staggered magnetic 
field (0.001 < h < 0.01), the values of the energy gap 
coming from a finite system basically represent the per- 
turbative behavior. We have shown that in this range of 
the staggered magnetic field, we are not allowed to read 
the scaling exponent of the energy gap directly from the 
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extrapolated numerical results. 

We have applied a general finite size scaling procedure 
for investigating the 7J-dependence of the critical expo- 
nent of the gap. We have identified two regimes through 
which the real critical exponent v{H) can be numerically 
calculated. To find the correct exponent of the gap in 
small- a; regime (x — Nh u ^ H ^ <C 1), we have used the 
scaling behavior of the coefficient of the leading term in 
the perturbation expansion, which is introduced by au- 
thors in Ref. [25] . In the large- x regime using the standard 
finite size scaling Eg. (fill)) , we have computed the correct 
critical exponent. On the other hand, using the confor- 
mal estimates of the small perturbation (h <C 1), we have 



found the independence of the critical exponent (Eq.©) 
from the theoretical point of view. Our numerical results 
in both regimes are in well agreement with the results 
obtained by the theoretical and numerical approaches. 



VI. ACKNOWLEDGMENTS 

I would like to thank A. Langari, M. R. H. Khajeh- 
pour, G. I. Japaridze, J. Abouie and M. Kohandel for 
insightful comments and fruitful discussions that led to 
an improvement of this work. 



* Electronic address: mahdavifarOiasbs.ac.ir 

1 D. C. Dender, P. R. Hammar, D. H. Reich, C. Broholm 
and G. Aeppli, Phys. Rev. Lett. 79, 1750 (1997). 

2 T. Asano, H. Nojiri, Y. Inagaki, J. P. Boucher, T. Sakon, 
Y. Ajiro and M. Motokawa, Phys. Rev. Lett. 84, 5880 
(2000). 

3 T. Asano, H. Nojiri, W. Hiegemoto, A. Koda, R. Kadono 
and Y. Ajiro, J. Phys. Soc. Jpn. 72, 594 (2002). 

4 M. Oshikawa, K. Ueda, H. Aoki, A. Ochiai and M. Kohgi, 
J. Phys. Soc. Jpn. 68, 3181 (1999). 

5 H. Shiba, K. Ueda and O. Sakai, J. Phys. Soc. Jpn. 69, 
1493 (2000). 

6 M. Kohgi, K. Iwasa, J. M. Mignot, B. Fak, P. Gegenwart, 
M. Lang, A. Ochiai, H. Aoki and T. Suzuki, Phys. Rev. 
Lett. 86, 2439 (2001). 

7 M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 
(1997); Phys. Rev. B 60, 1038 (1999). 

8 Y. -J. Wang, F. H. L. Essler, M. Fabrizio and A. A. Ners- 
esyan, Phys. Rev. B 66, 024412 (2002). 

9 M. Sato and M. Oshikawa, Phys. Rev. B 69, 054406 (2004) 

10 F. C. Alcaraz and A. L. Malvazzi, J. Phys. A 28, 1521 
(1995). 

11 R. Feyerherm, S. Abens, D. Gunther, T. Ishida, M. Meib- 
ner, M. Meschke, T. Nogami and M. Steiner, J. Phys. Con- 
dens. Matter 12, 8495 (200). 

12 V. L. Pokrovsky and A. L. Talapov, Soc. Phys. JETP 51, 
134 (1980). 



13 R. B. Griffiths, Phys. Rev. 133, A768 (1964). 

14 J. Lou, C. Chen, J. Zhao, X. Wang, T. Xiang, Z. Su and 
L. Yu, Phys. Rev. Lett. 94, 217207 (2005). 

15 A. Fledderjohann, M. Karbach and K.-H. Mutter, Eur. 
Phys. J. B 7, 225 (1999). 

16 A. L. Chernyshev, Phys. Rev. B 72, 174414 (2005). 

17 N. M. Bogoliubov, A. G. Izergin and V. E. Korepin, Nucl. 
Phys. B 275, 687 (1986). 

18 G. Grosso and L. Martinelli, Phys. Rev. B 51, 13033 
(1995). 

19 A. Fledderjohann, M. Karbach and K.-H. Mutter, Eur. 
Phys. J. B 5, 487 (1998). 

20 A. Klumper, Eur. Phys. J. B 5, 677 (1998), and references 
therein. 

21 F. C. Alcaraz, M. N. Barber, M. T. Batchelor, Phys. Rev. 
Lett. 58, 771 (1987). 

22 F. Woynarovich, H. -P. Eckle, J. Phys. A : Math. Gen. 20, 
L97 (1987). 

23 F. C. Alcaraz, M. N. Barber, M. T. Batchelor, Ann. Phys. 
182, 280 (1988). 

24 D. V. Dmitriev, V. Ya. Krivnov, A. A. Ovchinikov, and A. 
Langari, JETP 95, 538 (2002). 

25 A. Langari and S. Mahdavifar, Phys. Rev. B 73, 54410 
(2006). 

26 J. -B. Fouet, O. Tchernyshyov, and F. Mila, Phys. Rev. B 
70, 174427 (2004). 



